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ABSTRACT 

We present the results of a new search for variable stars in the Local Group (LG) isolated dwarf 
galaxy IC 1613, based on 24 orbits of F475W and F814W photometry from the Advanced Camera for 
Surveys onboard the Hubble Space Telescope. We detected 259 candidate variables in this field, of 
which only 13 (all of them bright Cepheids) were previously known. Out of the confirmed variables, 
we found 90 RR Lyrae stars, 49 classical Cepheids (including 36 new discoveries), and 38 eclipsing 
binary stars for which we could determine a period. The RR Lyrae include 61 fundamental (RRa6) 
and 24 first-overtone (RRc) pulsators, and 5 pulsating in both modes simultaneously (RRci). As for 
the majority of LG dwarfs, the mean periods of the RRab and RRc (0.611 and 0.334 day, respec- 
tively) as well as the fraction of overtone pulsators (f c =0.28) place this galaxy in the intermediate 
regime between the Oosterhoff types. From their position on the period-luminosity diagram and 
light-curve morphology, we can unambiguously classify 25 and 14 Cepheids as fundamental and first- 
overtone mode pulsators, respectively. Another two are clearly second-overtone Cepheids, the first 
ones to be discovered beyond the Magellanic Clouds. Among the remaining candidate variables, five 
were classified as 5-Scuti and five as long-period variables. Most of the others are located on the 
main-sequence, the majority of them likely eclipsing binary systems, although some present variations 
similar to pulsating stars. We estimate the distance to IC 1613 using various methods based on the 
photometric and pulsational properties of the Cepheids and RR Lyrae stars. The values we find are 
in very good agreement with each other and with previous estimates based on independent meth- 
ods. When corrected to a common reddening of E(B— V)=0.025 and true LMC distance modulus of 
(m— M) £mco=18.515±0.085, we find that all the distance determinations from the literature converge 
to a common value of (m-M) =24.400±0.014 (statistical), or 760 kpc. The parallel WFPC2 field, 
which lies within three core radii, was also searched for variable stars. We discovered nine RR Lyrae 
stars (4 RRa&, 4 RRc and l RRci) and two Cepheids, even though the lower signal-to-noise ratio of 
the observations did not allow us to measure their periods as accurately as for the variables in the 
ACS field-of-view. We provide their coordinates and approximate properties for completeness. 
Subject headings: binaries: eclipsing — galaxies: dwarf — galaxies: individual (IC 1613) — Local 
Group — stars: variables: Cepheids — stars: variables: RR Lyrae 



1. INTRODUCTION 

The search for variable stars in IC 1613 started at the 
beginning of the 1930's with the seminal work of Baade, 
H ubble, and May all, although it was published only later 
m ISandagd (fl97llh The discovery of 37 definite Cepheids 
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allowed Baade to calculate the first distance to IC 1613 
based on the period-luminosity (PL) relation. The analy- 
sis of the faint Cep heids in the same pho t ograp hic plates 
was completed by iCarlson fc Sandagi (| 19901) and the 
PL relation extended to periods as short as about two 
days. Candidate RR Ly rae stars were a lso observed 
in IC 1613 fro m ground- dSaha et al.l 119921 ) and space- 
based imaging (|Dolr>hin et al.ll 200D. confirming the pres- 
ence of a bona fide old popula tion. More recent sur- 
veys by IMantcgazz a et ahl (|2001|, and references th erein) 
and the OGLE collaboration (|Udalski et al.l 120011) of a 
central ~14'xl4' field led to the discovery of 128 and 
138 Cepheids, respectively, although the small size of 
the telescopes used by both groups limited the sample 
to relatively bright, long-period Cepheids. While all 
the Cepheids with periods longer than about 10 days 
in IC 1613 have most likely been discovered, very few 
have been found at the faint end of the PL relation (e.g., 
iDolphin et~aDl2001h . 

Here we present the analysis of very deep Hubble Space 
Telescope (HST) ACS data with the aim of searching for 
short-period variable stars. These data are part of a 
larger set obtained in the context of the LCID project 
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(C. Gallart et al. 2010, in preparation). The depth 
and quality of the data also allow the reconstruction of 
the full star formation history (SFH) of this field located 
just outside of the core radius. The study of the stellar 
populations and SFH from the color-magnitude diagram 
(CMD) fitting technique will be presented in a compan- 
ion paper (E. Skillman et al. 2010, in preparation). 

The structure of the p resent paper follows closely that 
of iBernard et al.l ((20091 hereafter Paper I) in which we 
analysed the properties of the variables in the dSph 
galaxies Cetus and Tucana. A summary of the obser- 
vations and data reduction is presented in § [21 while § [3] 
and § S] deal with the identification of variable stars and 
the sample completeness, respectively. In §[5] we describe 
the sample of RR Lyrae stars and the Cepheids are pre- 
sented in § IH1 The eclipsing binaries and other candidate 
variables are presented in §[7]to[5J Our list of variables is 
cross-identified with the previous catalogs covering this 
field in §[TUJ In §[TT]we use the properties of the RR Lyrae 
stars and Cepheids to estimate the distance to IC 1613 
and compare it with the values of the literature. In the 
last section we discuss the results and present our con- 
clusions. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Primary ACS Imaging 

This work is based on observations obtained with the 
ACS onboard the HST of a field located about 5' West 
of the center of IC 1613. The field location was chosen to 
be far enough out in radius to minimize crowding effects 
for high quality photometry, yet close enough in to max- 
imize the total number of stars observed. As the goal of 
these observations was to reach the oldest main sequence 
turn-offs with good photometric accuracy, 24 HST orbits 
were devoted to this galaxy. These were collected over 
about 2.4 consecutive days between 2006 August 28 and 
30. Each orbit was split in two ^1200 seconds exposures 
in F475W and F814W for an optimal sampling of the 
light curves. The complete observing log is presented in 
Table [U 

The DAOP HOT / ALLFRAME suite of programs 
(|Stetsonlll99l was used to obtain the instrumental pho- 
tometry of the stars on the individual, non-drizzled im- 
ages provided by the HST pipeline. Additionally, we used 
the pixel area maps and data quality masks to correct for 
the variations of the pixel areas on the sky around the 
field and to flag bad pixels. Standard ca li bratio n was 
carried out as described in ISirianni et al.l ( 20051). tak- 
ing in to account the updated zero-points of iMack et al.l 
(|2007f ) following the lowering of the Wide Field Channel 
temperature setpoin t in July 2006. We refer the reader to 
IMonelli et al.l (|2010T > for a detailed description of the data 
reduction and calibration. The final CMD is presented in 
Fig.[]](>/t panel), where the (F475W + F8UW)/2 ~ V 
filter combination was chosen for the ordinate axis so 
that the horizontal-branch (HB) appears approximately 
horizontal. 

The light-curves of the Cepheids and RR Lyrae stars 
were converted to Johnso n BVI m agnitudes using the 
transformations given in IPaper J . Since 12 of our 
Cepheids were also observed by the OGLE team in V and 
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Dwarfs: 



a Modified Heliocentric Julian Date of mid-exposure: HJD — 2,400,000. 

I (see § [TO]) , we could check the accuracy of our trans- 
formations. Ten out of the 12 stars in common have 
well sampled light-curves and accurately determined pe- 
riods, although 6 of these only have V observations in 
the OGLE database because of their relatively low lumi- 
nosities. 

We find excellent agreement at all phases between the 
two photometry sets for 7 of the 10 usable stars. The re- 
maining three present a shift in magnitude which ranges 
from ~0.1 at maximum light to ~0.5 for the faintest 
points. The comparison of the finding charts shows that 
this difference is due to contamination by nearby bright 
stars that are not resolved in the ground-based data. To 
illustrate this, in Fig. [2] we compare the finding-charts 
and photometry of the OGLE database with ours for 
two Cepheids (V106 and V154). For each Cepheid, we 
show 15"xl5" stamps from the OGLE I-band and our 
F814W images, as well as the V and I light-curves from 
both photometries. While both Cepheids are located in 
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Fig. 1. — Color-magnitude diagrams of IC 1613 for the ACS (left) and WFPC2 (right) fields, where th e confirmed and candida te variables 
have been overplotted, as labeled in the inset. An isochrone from the BaSTI library (Z=0.0003, 13Gyr, Pictrinferni et al. 2004) is overlaid 
in the right panel to show the location of the old main-sequence turn-off. 



rather crowded fields, the stars close to V154 are rela- 
tively faint compared to the Cepheid. The ground-based 
and HST light-curves are therefore overlapping perfectly. 
On the other hand, in the case of V106 the neighbours 
have a luminosity similar to that of the Cepheid, thus af- 
fecting the mean magnitude and amplitude of the OGLE 
light-curve. 

This test shows that our photometry in the Johnson 
bands, although transformed from a very different pho- 
tometric system, is very reliable and can be safely used to 
analyse the properties of the variable stars, both in terms 
of magnitude and amplitude. This is emphasized in § 1111 
where the distance to IC 1613 is calculated from Johnson 
luminosities and for which we find excellent agreement 
with previous measurements based on other methods. 

2.2. Parallel WFPC2 Imaging 

IC 1613 was also observed with the WFPC2 in the 
F450W and F814W bands as parallel exposures to the 
primary ACS observations. This provided the opportu- 
nity to analyze a second field with the same number of 



TABLE 2 

Photometry of the Variable and Candidate Variable 
Stars in IC 1613 - HST Bands 



MHJD a 
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0.052 



Note. — Tabic f2] is published in its entirety in the electronic edi- 
tion of The Astrophysical Journal. A portion is shown here for guid- 
ance regarding its form and content. 



a Modified Heliocentric Julian Date of mid-exposure: 
H,TD-2,400,u00. 

observations and a similar exposure time. The orienta- 
tion of the ACS field was chosen such that the parallel 
WFPC2 field would sample a region further from the 
center of the galaxy (~11') in order to study population 
gradients. The images from the Wi de Field chips were 
reduced individually as described in [Turner] (1997), and 
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TABLE 3 

Photometry of the Variable and Candidate Variable Stars in IC 1613 - Johnson Bands 



MHJD a 


B 




MHJD a 


V 
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07 
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24.654 


0.052 



Note. — Tablc[3]is published in its entirety in the electronic edition of The Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 



a Modified Heliocentric Julian Date of mid-exposure: HJD — 2,400,000. 



TABLE 4 

Properties of Variable Stars in IC 1613 



ACS Field 







R.A. 


Decl. 


Period 
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Note. 



Table [4*1 is published in its entirety in the electronic edition of The Astrophysical Journal. A portion is shown here for guidance regarding its form a 



calibrated following the instructions from the ffiSTData 
Handbook for WFPC2F ro l Given the small number of 
stars on the Planetary Camera, the photometry of this 
chip was performed on the stacked images only. The re- 
sulting CMD for the whole WFPC2 field-of-view is shown 
in the right panel of Fig. [H where a 13 Gyr old isoc hrone 
from the BaSTI library (jPietrinferni et all I2004D has 
been overplotted as suming Ab=0.123 and A/=0.056 
(|Schlegel et alJll998ft for the F450W and F814W bands, 
respe ctively , and a dereddened distance modulus of 24.40 
(see S HOt . 

Unfortunately, due to the lower sensitivity of the in- 
strument with respect to the ACS, especially in the 
F450W band, the CMD does not reach the oldest main- 
sequence turn-offs. However, the lack of bright main- 
sequence stars in this outer field reveals the presence of 
a blue horizontal-branch component, typical of old, low 
metallicity stellar populations. 

3. IDENTIFICATION OF VARIABLE STARS 

The candidate variables were extracted from the 
photometric cat a logs using the variability index of 
iWelch fe Stetson! (jl993f ). This process yielded ~780 can- 
didates in the primary field. A preliminary check of the 
light-curve and position on the CMD, together with a 
careful inspection of the stacked image, allowed us to dis- 
card false detections due to cosmic-ray hits, chip defects 
or stars located under the wings of bright stars. This 
left us with 259 candidates; these are shown in Fig. [1] 
overplotted on the CMD using their intensity-averaged 
magnitudes when available (see below). The individual 
F475W and F814W measurements for all of the candi- 
date variables are listed in Table [2j while the transformed 
BVI magnitudes are given in Table [3] 

The period search was performed on the suspected 
variables through Fourier analysis (jHorne fc Baliunasl 

10 http:/ /www. stsci.edu/instruments/wfpc2/Wfpc2_dhb/wfpc2_ 
ch52.html 



1986) taking into account the information from both 
bands simultaneously, then refined interactively by mod- 
ifying the period until a tighter light curve was obtained. 
For each variable, datapoints with error bars larger than 
3-<7 above the mean error bar size were rejected through 
sigma clipping with five iterations. As the period-finding 
program is interactive, it was possible to selectively re- 
ject more or less datapoints depending on the light curve 
quality before recalculating the periodogram. However, 
except in few particular cases, we found that the period- 
search was not affected by a few bad points. Given the 
short timebase of the observations (<3 days), the pe- 
riods are given with three significant figures only. The 
accuracy of period determination was estimated from the 
degradation of the quality of the light curves when ap- 
plying small period offsets. It mainly depends on the 
period itself and on the time interval covered by obser- 
vations, and ranges from about 0.001 day for the shorter 
period RR Lyrae stars to few hundredths of a day for the 
longest period Cepheids. 

The classification of the candidates was based on their 
light-curve morphology and position in the CMD. We 
found 90 RR Lyrae stars, 49 Cepheids, 38 eclipsing bi- 
naries, and five candidate 5-Scuti stars. A significant 
number of variable candidates was also found along the 
main sequence. For most of these, however, we could not 
find a period that would produce a smooth light-curve. 
In the following, we will simply refer to these variables as 
main-sequences variables (MSV), although most of them 
are probably eclipsing binaries. 

To obtain the amplitudes and intensity-averaged mag- 
nitudes of the monoperiodic variables in the instability 
strip, we fitted the l ight-curves wi t h a se t of template s 
based on the set of lLavden et al.1 (|1999t see iPaper ID . 
The amplitudes of the double-mode RR Lyrae stars were 
measured from a low-order Fourier fit to the light-curve 
phased with the primary period after prewhitening of 
the secondary period. The mean magnitude and color 
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TABLE 5 

Candidate Variable Stars in IC 1613 - ACS Field 







R.A. 


Decl. 






ID 


Type 


(J2000) 


(J2000) 


F475W 


F814W 


VC194 


MSV 


01 04 20.27 


02 10 28.6 


25.272 


25.261 


VC195 


RGV 


01 04 20.62 


02 10 01.1 


25.777 


24.361 


VC196 


MSV 


01 04 22.63 


02 09 53.7 


26.401 


26.078 


VC197 


MSV 


01 04 23.08 


02 08 06.6 


26.680 


25.784 


VC198 


MSV 


01 04 23.24 


02 08 33.6 


24.092 


24.221 



Note. — Tabic ]E\ is published in its entirety in the electronic 
edition of The Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 



of the variables for which we could not find a period 
are weighted averages from the ALLFRAME photome- 
try and are therefore only approximate. Table 2] sum- 
marizes the properties of the confirmed variable stars in 
the ACS field. The first two columns give the identifica- 
tion number and type of variability, while the next two 
list the equatorial coordinates (J2000.0). Columns (5) 
and (6) give the primary period in days, i.e., the first- 
overtone period in the case of the RRd, and the loga- 
rithm of this period. The intensity-averaged magnitudes 
(FA75W) and (F8UW), and color (F475W) - {F8UW) 
are given in columns (7), (9), and (11). The amplitudes 
in the F475W and F814W bands measured from the tem- 
plate fits are listed in the eighth and tenth columns. The 
last six columns alternately list the intensity-averaged 
magnitudes and amplitudes in the Johnson B, V, and I 
bands. Approximate values are listed in italics. A list 
of the remaining candidates, providing the coordinates 
and approximate magnitudes, is given in Table [5j They 
are labeled as MSV, long-period variables (LPV), classi- 
cal instability strip variables (ISV), and RGB variables 
(RGV) based on their location on the CMD. 

The same procedure was followed with the WFPC2 
photometry, in which we found 11 variables: nine 
RR Lyrae stars and two Cepheids. However, the low 
signal-to-noise of the individual measurements produced 
rather noisy light-curves, and some low amplitude vari- 
ables might have been missed. Given the small number 
of variables in the parallel field and the generally lower 
quality of their photometry and inferred parameters, we 
only give their coordinates and approximate parameters 
in Table |6] for completeness. These variables will not be 
taken into account when calculating average properties 
of the IC 1613 variable star population. 

4. COMPLETENESS AND AREAL COVERAGE 

In this section, we analyze the effects of stellar crowd- 
ing and signal-to-noise (SNR) limitations, temporal sam- 
pling, and spatial coverage on the completeness of our 
Cepheids and RR Lyrae stars samples. 

The high spatial resolution of the ACS and the depth of 
our data imply that incompleteness will become notice- 
able well below the HB. Artificial-star tests (E. Skillman 
et al. 2010, in preparation) indicate that the complete- 
ness is higher than 98% at (F475W+F814W)/2 - 25.0. 
Therefore, down to these magnitudes only variables with 
amplitude of the order of the error bars at this mag- 
nitude (~0.1) might have been missed. However, vari- 
ables fainter than the HB (e.g., SX Phoenicis and/or 
(5-Scuti) have probably been missed due to crowding and 
low SNR. In addition, even though these variables are 




V154 (OGLE 2760) 



0.0 0.5 1.0 

Phase 




V106 (OGLE 21 1 7) 

23.5 1 , , i , , , , i , , : 

0.0 0.5 1.0 

Phase 

Fig. 2. — Comparison of the finding charts and Johnson photom- 
etry from the OGLE database with ours for V154 {top) and V106 
{bottom). The finding charts arc 15" on a side. The OGLE light 
curves are represented as open symbols, while our photometry is 
shown as filled symbols. For clarity, the V-band datapoints have 
been shifted downward by 0.2 mag. Note the good agreement of 
the ground-based and /^ST light-curve in the first case. The second 
is an example of ground-based photometry affected by unresolved 
bright companions, which affects the light-curve morphology and 
mean- magnitude. 

present in this galaxy (see § [8]) , the relatively long expo- 
sure time smoothing out the variations in luminosity and 
the rather slow temporal sampling hampered the detec- 
tion of these short-period variables. 

We also estimate the completeness due to temporal 
sampling by carrying out numeric al simulations in the 
same way as described in iPaper 1 The detection prob- 
ability is presented in Fig. [3] In the upper panel, the 
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TABLE 6 

Properties of Variable Stars in IC 1613 - WFPC2 Field 







R.A. 


Decl. 


Period 








(F450W0- 






ED 


Type 


(J2000) 


(J2000) 


(days) 


log P 


(F450W) 


<F814iy) 


{F81AW) 


^4-450 


^4-814 


V183 


d 


01 04 04.90 


02 13 44.9 


0.34 


-0.46 


25.178 


24.446 


0.732 


0.463 


0.146 


V184 


c 


01 04 05.05 


02 13 11.0 


0.36 


-0.44 


24.955 


24.313 


0.642 


0.520 


0.299 


V185 


ab 


01 04 05.20 


02 13 32.7 


0.56 


-0.25 


25.150 


24.240 


0.910 


1.168 


0.551 


V186 


c 


01 04 05.29 


02 11 01.7 


0.35 


-0.46 


25.107 


24.272 


0.834 


0.510 


0.234 


V187 


ab 


01 04 07.61 


02 13 00.8 


0.55 


-0.26 


25.096 


24.393 


0.703 


1.615 


0.565 


V188 


ab 


01 04 07.99 


02 11 32.5 


0.65 


-0.19 


25.097 


24.174 


0.924 


0.794 


0.409 


V189 


ab 


01 04 08.65 


02 11 11.0 


0.60 


-0.22 


25.228 


24.279 


0.949 


0.773 


0.296 


V190 


Cep 


01 04 10.25 


02 13 06.9 


0.67 


-0.17 


23.507 


22.693 


0.813 


0.664 


0.401 


V191 


c 


01 04 10.32 


02 12 50.9 


0.33 


-0.48 


24.979 


24.359 


0.619 


0.624 


0.344 


V192 


Cep 


01 04 10.40 


02 12 03.7 


1.30 


0.11 


22.701 


21.896 


0.805 


0.657 


0.285 


V193 


c 


01 04 10.97 


02 12 08.1 


0.36 


-0.44 


25.072 


24.411 


0.661 


0.516 


0.276 
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Fig. 3. — Probability of detecting variable stars in IC 1613 as 
a function of period, for periods between about 1 hr and 5 days 
(top), and close-up view for periods between about 1 hr and 1 day 
(bottom). 



maximum period that is displayed is limited by the ob- 
servational timebase (~3 days). Note that even though 
it is possible to detect variable stars with periods longer 
than the observation time-base, it is not feasible to de- 
termine their periods. The Figure also shows that the 
probability to detect a variable star in the period range 
of RR Lyrae stars (~ 0.2-0.8 day) is basically one. On 
the other hand, a significant dip is present at P~l day, 
due to the observations being gathered at the same time 
of the day each day. While it might not prevent detect- 
ing variables at this period (most likely Cepheids) if their 
amplitude is larger than about ^0.1, it seriously limits 
our capacity to find an accurate period (e.g., V150 in 
Fig.H V030 in §iU). 

As expected from the large extent of the galaxy com- 
pared to the field-of-view of the ACS, the main factor pre- 
venting us from obtaining a complete census of RR Lyrae 
stars and Cepheids in IC 1613 is the spatial coverage. 
Figure 0] shows the distribution of stars from our ACS 
and WFPC2 fields on to p of ellipses representi ng the core 
radius (r c = 4.5 ' ± 0.9. IBattinelli et~aT1l2007l) and twice 
the core radius. IBattinelli et al.l (|2007|) estimated a tidal 



radius of 24.4' ±0.3, ellipticity e = 1 -b/a = 0.19 ±0.02, 
and position angle of 87°, in g ood a greement with the val- 
ues found bv lBernard etall ([20071 PA = 80°, e = 0.15). 
Thus, the ACS field covers only about l/160th of the 
area within the tidal radius. If we assume that the dis- 
tribution of variable stars follows that of the overall pop- 
ulation, we can estimate the total number of RR Lyrae 
st ars and Cepheids within the tidal radius as described 
in lPaper H Basically, we adopt the shape and orientat ion 
of the isodensity contours from IBattinelli et al.l (|2007[ ) to 
divide our sample of variable stars into elliptical annuli. 
The area of these annuli was obtained through Monte- 
Carlo sampling, from which we could calculate the den- 
sity profile of the variable stars. Using thi s profi le with 
the core and tidal radii of IBattinelli et ail ([20 07) as in- 
put to equation (21) of iKingl (j!962l ). we estimate that 
there should be about 2100 RR Lyrae stars and 1400 
Cepheids within the tidal radius of IC 1613. Therefore, 
our ACS sample represents about 4% of the total number 
of RR Lyrae stars and Cepheids. 

5. RR LYRAE STARS 
5.1. ACS sample 

In our ACS field we detected 90 RR Lyrae stars, all 
of them new discoveries (see § I5.3j) . From the periods 
and light-curve shapes, we identified 61 of them pulsat- 
ing in the fundamental mode (RRa6), 24 in the first- 
overtone mode (RRc), and five in both modes simultane- 
ously (RRd). The light curves of all the RR Lyrae stars 
are shown in Fig. [5] 

We find the mean period of the RRa6 and RRc to be 
0.611 and 0.334 day, respectively. These values, as well as 
the fraction of first-overtone f c = N c / (N a b + N c ) = 0.28, 
place IC 1613 in the gap between the Oosterhoff type I 
and type II groups (Oo-I and Oo-II). IC 1613 there- 
fore joins the majority of dwarf galaxies in the Local 
Group as an Oosterhoff intermediate galaxy, the only 
real exceptions among th e dSphs with well-s ampled HB 
to date being Ursa Mi nor (jNemec et al.lll988() and Bootes 
(IDairOra et al.ll200l . 

In Fig. [6] we show the period- amplitude diagram (top) 
and the period distribution (bottom) of the RR Lyrae 
stars. The top panel shows that the distribution of the 
RRafr of IC 1613 in period-amplitude space is very similar 
to that of Tucana, both in terms of slope and dispersion 
around the fit. The slight shift towards shorter periods of 
the RRa& stars of IC 1613, compared to those of Tucana, 
might be explained by the higher metallicity of IC 1613. 
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Fig. 4. — Spatial distribution of stars in the ACS and WFPC2 fields, wher e the variable stars ha ve been overplotted as labeled in the 
inset. Solid- and dashed-line ellipses represent the core radius (r c = 4.5' ± 0.9, Battinclli ct al. 2007) and twice the core radius. 



However, the distribution of fundamentalized periods 
appears to be bimodal, as shown by the fit of two Gaus- 
sian profiles (solid grey line) in the bottom panel, with 
peaks at P=0.478 and 0.614 day. Given the separa- 
tion between the two peaks and the small uncertainty 
on the periods (<0.001 day), it seems unlikely to be 
caused by stochastic effects only. To check the relia- 
bility of the observed bimo dality, we applied the KMM 
test of lAshman et al.l ([1994) to the sample of RRafe peri- 
ods. The algorithm estimates the improvement of a two- 
group fit over a single Gaussian, and returned a p-value 
of 0.075. This can be interpreted as the rejection of the 
single Gaussian model at a confidence level of 92.5%, and 
therefore indicates that a bimodal distribution is a statis- 
tically significant improvement over the single Gaussian. 

Note that the distribution of periods of the RRafr alone 
is also bimodal, although the small number of RRafr 
in the short-period peak does not exclude a result due 
to small number statistics. In addition, the location of 
the peaks (at P=0.494 and 0.615 day) does not seem to 
correspond to the superposition of Oo-I and Oo-II sam- 
ples, as is the case in NGC 1835 (P~0.54 and 0.65 day; 
iSoszvriski et~aT1[200l . 

The bimodality might nontheless be due to the pres- 
ence of two distinct old populat ions, as was alread y 
observed in the case of Tucana ([Bernard et al.l 120081) . 



Figure [7] presents a zoom-in of the CMD centered on 
the HB. Zero-age horizontal-branch (ZAHB) loci and 
tracks from the scaled-s olar models of the BaSTI library 
(jPietrinferni et aLl 120041 ) are also shown, plotted assum- 
ing a reddening of E(B— V)=0.025 from ISchlegel et all 
(1998), and a distance modulus of 24.49. This dis- 
tance takes into account the true distance modulus 
(m— M)o=24.44 determined in ij |11.3l based on these data 
and a shift of +0.05 t o correct for the up dated electron- 
conduction opacities (jCassisi et al.l 12007). This Figure 
shows that more than half of the RR Lyrae stars are con- 
centrated in the faintest ~0.1 mag of the instability strip 
(IS), while the other half is spread over ~0.3 mag above. 
From the evolutionary tracks, it appears that most of the 
variables belong to the fainter, more-metal rich sample; 
these stars tend to have shorter periods. The brighter 
variables, on the other hand, appear to be related ei- 
ther to the Z=0.0001 population, or to stars with slightly 
higher metallicity (Z~0.0003) that have evolved off the 
blue horizontal-branch visible on the WFPC2 CMD in 
Fig.ffl 

In any case, the range of luminosity covered by the HB 
stars within the IS cannot be explained by a monometal- 
lic population, even taking into account evolution off 
the ZAHB. Given that RR Lyrae stars are old stars 
(>10 Gyr), this implies a relatively quick chemical en- 
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Fig. 5. — Light-curves of the RR Lyrae variables of the ACS field in the F475W (black) and F814W (grey) bands, phased with the period 
in days shown in the lower right corner of each panel. Photometric error bars are shown. The open circles show bad data points, i.e., with 
errors larger than 3-<r above the mean error of a given star, which were not used in the calculation of the period and mean magnitudes. 
[FiaureUJlis presented in its entirety in the electronic edition of the Astrophysical Journal]. 



richment in the first few billion years of the formation of 
the galaxy. 

5.2. WFPC2 sample 

In the parallel WFPC2 field we discovered nine RR 
Lyrae stars (4 RRa6, 4 RRc and 1 RRcQ, even though 
the low signal-to-noise ratio of the observations did not 
allow us to measure their periods as accurately as for the 
variables in the ACS field-of-view. Based on the small 
samples of RR Lyrae stars, we find (P a b)=0.59±0.03 and 
(P c }=0.35±0.01 for the fundamental and first-overtone 



pulsators, respectively, which is in rough agreement with 
the values found for the RR Lyrae stars in the ACS field- 
of-view. We provide their coordinates and approximate 
properties for completeness in Table [51 

5.3. Comparison with known RR Lyrae stars in IC 1613 

Prior to our study, two surveys for variable stars found 
candidate RR Lyrae stars in IC 1613, although none 
of them in common with those presented here as the 
pointings did not overlap. U s ing th e 4-shooter on the 
Hale-5m telescope. iSaha et al.l(|1992l ) found 15 candidate 
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Fig. 6. — Top: Period-amplitude diagram for the RR Lyrae stars. 
Circles, triangles and squares represent RRab, RRc, and RRd (plot- 
ted with their overtone periods) respectively. The dashed line is a 
fit to the period-amplitude of the RRab after rejecting the points 
further than 1.5-a (dotted lines). The light and dark greyed ar- 
eas r epresent the ±1.5-cr limits of Cetus and Tucana, respectively, 
from IPaper J . Bottom: Period histogram for the RR Lyrae stars 
of the ACS field. RRab and RRc are shown as the light and dark 
gray histograms, respectively, while the thick black histogram rep- 
resent the fundamentalized RR Lyrae stars (RRab, RRc, and RRd) . 
The solid gray line is a double-Gaussian fit to the fundamentalized 
histogram. 

RR Lyrae stars in a 8'x8' field of view located 12' to the 
west of the galaxy center. From their Table 3, we calcu- 
lated a mean period for the RRa& of 0.60±0.02 day, which 
is in good agreement with our value. We note, however, 
that all their RR Lyrae stars are of the ab type and on 
average appear brighter by about 0.3 mag than expected 
from the distance and metallicity of IC 1613 (see also 
§ ITT]) . While some RRc stars could have been missed 
due to their low amplitude and the relatively noisy light- 
curves, some of them might have been misidentified as 
RRafr because of the small number of datapoints. Alter- 
natively, given that crowding seems to be insignificant 
in this field — according to the finding charts — and that 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

F475W-F814W 

Fig. 7. — Zoom-in on the HB of IC 1613, where the RR Lyrae vari- 
ables have been overplotted. Open circles, open triangles, and filled 
squares represent RRab, RRc, and RRd, respectively. The long- 
dashed, thick lines represent, from top to bottom, the Z=0.0001, 
0.0003, and 0.001 ZAHB, while the light gray, dark gray, and 
black lines are the corresponding evolutionary tracks with masses 
as labeled. Overplotted on them, the crosses indicate intervals of 
10 Myr. The vertical dashed lines roughly delimit the instability 
strip. 

their RRa& span a range of luminosity of over a magni- 
tude, we suspect that some of their RR Lyrae candidates 
could actually be short-period Cep heids. 

The analysis of WFPC2 data bv lDolphin etaLl IpOOlf l 
in a field 7.4' southwest from the center, on the other 
hand, returned 13 RR Lyrae candidates. Assuming that 
the likely and possible overtone pulsators are actual RRc, 
their field contains 9 RRa6 and 4 RRc with mean periods 
of -0.58±0.02 and -0.37±0.03 day, respectively. This 
is very similar to the values we found for our WFPC2 
field. The slightly shorter mean period for the RRafe, 
when compared to our ACS sample, may be due to the 
difficulty of detecting the low amplitude, longer period 
variables because of the lower quality of the data. 

In any case, these samples at larger radii are not large 
enough for a robust comparison with the ACS sample in 
order to reveal potential gradients in the properties of 
the RR Lyrae s tars a s was done for the Tucana dwarf in 
Bernard etall (|2008D . In the case of Tucana, the ACS 
field of view covered a larger fraction of the radius of the 
galaxy, allowing a radial gradient study within the ACS 
field alone. 

6. CEPHEIDS 

In the search for variable stars presented here we found 
44 variables with properties of Cepheids: about 2 mag 
brighter than the HB, in or close to the IS, and P>0.6 
days. Their light-curves are shown in Fig. |8l Another 
five were found for which we could not determine their 
period and/or intensity-averaged mean magnitude due 
to inadequate temporal sampling. These are discussed 
in § EH 

Figure [9] shows the period— W^u diagram for the 
Cepheids in IC 1613. Because Cepheids from the dif- 
ferent modes of pulsation have a different magnitude at 
a given period, their location in the period-luminosity 
(PL) diagrams can be used to differentiate their type. 
These PL relations are fundamental in that they repre- 
sent the most robust method to calculate the distance to 



10 



Bernard et al. 



21 
22 
23 
24 
21 
22 
23 
24 
21 
22 
23 
24 
21 
22 
23 
24 
21 
22 
23 
24 



i V006 p=0.680 


••* ••• 

V008 p=1.690 


'\ 

..... •• 

" ••••• 

V011 p=1.487 


— i — e — j — i i — ' ' — h- 1 — ' ' — i — 3 — 

• • ' 



• • '"*•• 

' V022 p=1 .480 " 


— i — i — | — i i — . . — E — i — > • — > — 3 — 

V033 p=1 .480 " 


— i — i — | — i i — . . — E — i — > ■ — i — 3 — 

• *•• • m 

# ^ ■ # ^ 

V035 p=1.055 ] 


— I 1 1 1 1 — 1 1 — | — H 1 ' — >H — ' ' — 

- - 

' V037 p=1 .440 " 


— ) — i — | — i i — i i — | — i — i i — H — >— < — 

■■-*-..w...-. " 
" V041 p=1.240 ~ 


— I — i — | — ' ' — ' ' — 1 — i — ' ' — >H — 

*• •• 
>«,.. . . 

V043 p=1 .600 } 


— < — —3 — ' ' — > • — — > ■ — i — 3 — 

. .■*"••• 

'•• „ * * • .. , ! 

. • * * 1 . • * " 1 

■• "(•■•■•■a* * 

" V054 p=0.633 " 


— I 1 1 !— I 1 i E— ( E— , 1 1 E— 

• ". • *. 



V057 p=1.690 " 


— < — ( — j — t i — ' ' — E — i — •— ■ — i — 3 — 

_ *".-....." 

»*•»' ' •*--.».*.** 

' V061 p=1.245 \ 


— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — 
•••••.„• 

*"*••• $ *"••••" 

: V068 p=1 .900 " 


— ! 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

'"• " 

" V077 p=1 .350 " 


— i — i — i — i — i — i — i — i — i — i — i — i — i — t— 
a *.....••■* 

.'i 

•••• 

V078 p=1 .400 i 



0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 
Phase Phase Phase 



Fig. 8. — Light-curves of all the Cepheid variables in the F475W (black) and F814W (grey) bands. Photometric error bars are shown. 
The open circles show the data points with errors larger than 3-<r above the mean error of a given star. [Figure\S[is presented in its entirety 
in the electronic edition of the Astrophysical Journal]. 



galaxies within the nearby universe. To reduce the scat- 
ter due to interstellar extinction and the intrinsic disper- 
sion due to the finite width of the instability-strip, and 
to obtain a more secure classification, we calculated the 
Wesenheit, or r edden ing-free, magnitude introduced by 
Ivan den Berghl (|1968f ) . In the HST bands used here, the 
relation is of the form: 

W sli = F8UW - 0.945(.F475VF - F8UW) (1) 

where the coefficient 0.945 comes from the standard 
interstellar extinction curve dependence of the F814W 



magni tude on E(F475W-F814W), from ISchlegel et all 
(fl998h . 

Figure IH1 shows two almost parallel sequences, and cor- 
respond to the fundamental (F, open circles) and first- 
overtone (FO, filled circles) Cepheids, respectively. A few 
outliers are marked as open triangles (V006 and V091) 
and a fi lled square (V126), and are discussed in § 16.21 
and 16.31 respectively. From this plot, we find that 25 
Cepheids are pulsating in the fundamental mode, while 
16 fall on the first-overtone sequence. 

Figure [10] shows the PL diagram for the Johnson VI 
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Fig. 9. — Period- Wesenheit diagram for Ccphcids in IC 1613 (see 
text for details). The dash-dotted and dashed lines are linear fits 
to the fundamental and first-overton e Cepheids, re spec tively. T he 
individual Cepheids mentioned in § 16.11 16.21 16.31 and 111. 1.11 are 
labeled. 

bands obtained as desc ribed in iPaoer 1 as w ell as in the 
Wi band as denned in lUdalski et all (|1999cL Wj = I - 
1.55(V — /)). The new Cepheids of IC 1613 (open circles) 
are shown overplottcd on the classical Cepheids of the 
Large and Small Magellanic Clouds (LMC fc SMC) and 
of IC 1613 from the OGLE collaboration ([Udalski et al.1 
I2001L and references therein) as labeled in the inset. The 
apparent magnitudes were converted to absolute mag- 
nitude ass uming a distance modu lus of 18.515±0.085 to 
the LMC (Clcmcnt ini et al.ll2003|) . and a distance offset 
of 0.5 1 of the SMC relative to the LMC (jUdalski et al.l 
1999c). The Cepheids of IC1613 were shifted accord- 
ing to the distance determined in § QTJ It shows that 
the properties of the Cepheids discovered in our field 
are in excellent agreement with the properties of OGLE 
Cepheids, therefore confirming their classical Cepheid 
nature. 

However, one can see that the Cepheids in IC 1613 
(both ours and OGLE) are mainly located toward the 
short-period extremity of the fundamental and overtone 
sequences when compared to the Cepheids in the Mag- 
ellanic Clouds. Studies of the variable st ars in the very- 
meta l poor dwarf irreg ular galaxies Leo A (TDolphin et al.l 
12002ft and Sextans A (|Dolphin et al.ll2003ft also revealed 
large numbers of s hort period C epheids. This agrees with 
the suggestion by IBonol (|2003T ) that "the minimum mass 
that performs the blue loop [decreases with decreasing 
metallicity, which] means that metal-poor stellar systems 
such as IC 1613 should produce a substantial fraction of 
short-period classical Cepheids." 

Interestingly, the IC 1613 fundamental-mode Cepheids 
perfectly follow the various PL relationships defined by 
the Magellanic Clouds Cepheids over the whole range of 
periods, from ~1 to ~40 days (see Fig.QJJ. This suggests 
that the metallicity dependence of these relations, if any, 
must be very small at such low metallicities. 

6.1. The peculiar Cepheids V037 & V054 

Two of the Cepheids falling on the first-overtone se- 
quence in Fig. [9] present very peculiar light-curves. These 
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Fig. 10. — Period-luminosity diagrams for Cepheids in IC 1613. 
Our Cepheids are shown as open circles, while the symbols used 
for the OGLE Cepheids are as labeled in the inset. 

are shown again in Fig. 1111 plotted as a function of Julian 
date (left) and phase (right). 

In the top left panel, the amplitude of V054 seems to 
be decreasing at each consecutive maximum. The dis- 
persion at maximum light on the phased light-curves is 
reminiscent of the light-curves of KRd stars, although 
in this particular case the periodogram does not present 
the double-peak features characteristic of double-mode 
stars. However, the rather long main period of this vari- 
able (0.633 day vs. ~0.39 for KRd) limited the number of 
observed cycles to <4 and might explain the lack of res- 
olution in the periodogram. In addition, the position of 
V054 on the PL diagrams, at the shorter period end, is 
in excellent agreement with the double-mode Cepheids 
pulsating simultaneously in the fundament a l and first- 
overtone mode presented in ISosz vhski et al. (2008al). 

In the case of V037, we found that the only period 
giving a smooth light-curve is P=1.44 days. However, 
as shown in the bottom right panel of Fig. [Til this pro- 
duces a light-curve with a secondary maximum at phase 
~0.6. If this period is confirmed to be true, this would 
be the second Cepheid with such a light-curve in IC 1613 
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Fig. 11.— The peculiar Cepheids V054 (top) and V037 (bottom), 
plotted as a function of Julian date (left) and phase (right). 

after V39 discovered by Baade (see|Sandagc 197TJ) . Note, 
however, that more recent observ ations and analysis of 
V39 (e.g., Mantcga zza et ai1l200 If) favored a period only 
half that given by ISandagd 1197U P=14.350 vs. 28.72 
days) associated with a long term modulation of —1100 
days, which could correspond to an extremely distant 
W Vir star in the outer halo of the MW blended with a 
long-period variable in IC 1613. 

In our case, the short timebase of our observations def- 
initely rules out the possibility of long term modulation. 
In addition, a shorter period for V037 cannot phase the 
light-curve properly unless a strong modulation of the 
period itself is taken into account, which is very unlikely 
on such short timescales. Given that its position in the 
CMD and PL diagrams is exactly that expected from a 
Cepheid, in the following we will assume that V037 is a 
bona fide classical Cepheid. 

6.2. Second- Overtone candidates V006 & V091 

In the diagrams shown in Figs. |H] & EH two Cepheids 
lie above the first-overtone sequence in all panels. In 
particular in Fig. 1 1 01 they follow the sequence represent- 
ing the second-overtone pulsators found by the OGLE 
collaboration in the SMC. Their nearly sinusoidal light- 
curve with very low amplitude in both bands (<0.14 in 
F475W, -0.07 in F814W), as well as their location well 
within the IS, make them very strong second-overtone 
candidates. These are the first ones detected beyond the 
Magellanic Clouds. 

Second-overtone classical Cepheids are very rare ob- 
jects, as witnessed by the few firm candidates known 
to date: 1 3 out of 2049 classic al Cepheids — 0.42% — in 
the SMC (|Udalski et al.lll999aft. and 1 4 out of 3361— 
Q.63%— in the LMC (ISoszvnski et all l2008al). How- 
ever, theoretical models by lAntonello fc Kanburl ()1997f ) 
showed that these Cepheids are more likely to be found 
in low metallicity systems, which is in qualitative agree- 
ment with the fraction of second-overtone Cepheids in 
our IC 1613 sample, about 4%. 

6.3. V126: a blended Cepheid? 
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Fig. 12. — Sample light-curves for eclipsing binaries in IC 1613. 
The three brightest (top), three intermediate magnitude (middle), 
and the three faintest (bottom) binaries are shown. 

The light-curve morphology, luminosity and period of 
V126 are all characteristic of a classical Cepheid. On 
the other hand, its color is 0.5 mag bluer than the blue 
edge of the IS, and therefore appears below the funda- 
mental mode pulsator sequence in the period- Wesenheit 
diagram shown in Fig. [9] Although it appears isolated 
on the stacked images, the most likely hypothesis is that 
it is a bona fide Cepheid blended with an unresolved 
main-sequence star. The low amplitude of its luminosity 
variation in both bands relative to the other Cepheids 
with similar periods seems to support this conclusion. 

6.4. Other candidates 

Five Cepheids were observed for which we could not 
accurately determine their period and/or mean magni- 
tude due to inadequate temporal sampling. We use the 
weighted mean magnitude given by ALLFRAME as their 
approximate magnitudes. They are shown as open stars 
in Fig. [2 

V030 has a period —0.98 day, and was therefore ob- 
served at the same phase each day, at minimum light. 
Most of the rising phase and peak are missing so it was 
not possible to fit a template light-curve and measure 
accurately its mean magnitude. From its approximate 
magnitude, its location on the PL diagram is consistent 
with the fundamental mode Cepheids. 

V090 has relatively low-amplitude (A475— 0.15, 
Agi4~0.ll) and period 2<P<4 days. From its location 
on the PL diagram it is a probable first-overtone. 

VI 24 is very bright and was previously known from 
ground-based observations (see § [T0f . Given its very long 
period (—26 days), our observations cover less than 10% 
of a cycle so we did not try to obtain the mean magnitude 
from our data. It is shown in Fig. llOl as the fundamental 
mode OGLE Cepheid at log P=1.41. 

V137 was also present in the OGLE field, although 
their quoted period of 1.376 days does not provide a satis- 
factory fit to our data. Even though we could not obtain 



Faint Variables in IC 1613 



13 



an accurate period for this variable, the 1.4- mag decrease 
in the F475W band over the first two days and the begin- 
ning of the third, and the subsequent increase implies a 
period in the range 2.7<P<3.3 days, which would make 
it a fundamental mode Cepheid. 

VI 51 is one of the brightest Cepheids of our sam- 
ple, and thus has a period longer than our observa- 
tional timebase. Unfortunately, it is missing from the 
OGLE catalog so we can only constrain the period to 
the range 2.5<P<3.1 days from the light-curve morphol- 
ogy. It places it on the PL relation of fundamental mode 
Cepheids. 

7. ECLIPSING BINARY STARS 

A large number of variable star candidates were found 
in the main-sequence, as can be expected from the fact 
that it is well populated from the turn-offs to the B gi- 
ants. Most of the variables for which we could find a 
period are eclipsing binary stars. The remaining ones 
have properties of high-amplitude J-Scuti stars (HADS) 
and are discussed in the following section, while the can- 
didates for which the light-curves could not be phased 
are presented in § [9j Note that all of these stars are 
newly discovered variables. 

In the ACS field-of-view, we found 38 eclipsing bina- 
ries between 'V'~20.5 and ^26.5. At fainter magnitudes, 
the signal-to-noise ratio of the individual datapoints was 
too low to estimate the period and/or confirm variability. 
The light-curves of the three brightest, three intermedi- 
ate magnitude, and the three faintest binaries are shown 
in Fig. [H 

Given the small number of datapoints and relatively 
low SNR of most of the candidates, some periods might 
actually be multiples of the true periods. We thus did 
not try to assign them a particular type of binary (de- 
tached or contact) — although the variety of light-curve 
morphologies indicate that members of both kinds ap- 
pear to be present — nor check for their membership in 
IC1613. 

Bright eclipsing binaries have been used to determine 
distances to LG galaxies with claimed accurac y bet- 
ter than about 5% (e.g., M31: iRibas et all [20051 M33: 
iBonanos et al.l I2006T) . Unfortunately, the stars need to 
be bright enough (V ma x <20) to accomplish this in a 
reasonable time with current observing facilities. Even 
though some of our eclipsing binaries have accurate pe- 
riods and deep eclipses, none of them is brighter than 
V-20.4. 

8. CANDIDATE 5-SCUTI STARS 

Among the candidate variables fainter than the HB 
that appeared to be periodic variables, most of them 
were found to be binaries and are described in the pre- 
vious section. For some of these variables, however, it 
was possible to obtain a relatively smooth, pulsation-like 
light-curve (i.e., with a single minimum) with a short 
period. This was the case for 13 variables, which were 
flagged as possible J-Scuti stars. To further constrain 
their nature, additional parameters were taken into ac- 
count. First, we required that their color be within the 
classical IS, defined approximately by the position of the 
Cepheids and RR Lyrae stars. The top panel of Fig. Q2] 
shows a zoom-in of the CMD showing the location of the 
IS, where the 13 possible 5-Scuti are shown as filled sym- 
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Fig. 13. — Selection criteria for the 5-Scuti candidates. Top: 
Detail of the CMD, where the possible candidates are overplottcd 
with larger symbols. The dash-dotted lines indicate the approxi- 
mate location of the classical instability strip fitted on these data. 
Bottom: F814W versus F475W amplitudes, where the pulsating 
variables and binary stars have a different slope. Pluses, open cir- 
cles, and open stars are for the RRab, RRc, and Cepheids, respec- 
tively, while the open diamonds represent the binaries with good 
light-curves. The dotted and dash-dotted lines are linear fits to the 
pulsating variables and binaries, respectively. In both panels the 
5-Scuti candidates satisfying both requirements are shown as filled 
stars, and the rejected candidates as filled diamonds. 

bols. Eight of these fall within or close to the boundaries 
of the IS. 

Given that the pulsations of the 5-Scuti stars, like 
the Cepheids and RR Lyrae stars are driven by the re- 
mechanism, the second criteria was that the amplitude 
ratio had to be in agreement with that of the other pul- 
sating variables. The bottom panel of Fig. H31 shows the 
amplitude in F814W versus the amplitude in F475W for 
Cepheids, RR Lyrae stars, eclipsing binaries with good 
light-curves, and the possible 5-Scuti stars. Only 5 sat- 
isfy this requirement, and are shown as filled stars in both 
panels. They also present relatively smooth light-curves 
with very short periods (<0.2 days). Their pulsational 
properties are summarized in Table 0J 

9. OTHER VARIABLES 

In addition to the classical IS variables and eclipsing 
binaries presented above, we detected another 77 candi- 
date variables throughout the CMD, for which we could 
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Fig. 14. — Samp le li ght-curves for MSVs in IC 1613, as a function 
of HJD. As in Fig. 1121 the three brightest (top), three intermediate 
magnitude (middle), and the three faintest (bottom) MSV candi- 
dates are shown. Note the different vertical scale for the bottom 
panels. 

not determine the period because of inadequate tempo- 
ral sampling. Most of these candidate variables are lo- 
cated on the MS, although five LPV candidates are also 
present at or close to the TRGB. While the majority of 
the MS candidates are most likely eclipsing binaries, at 
least in the brightest part of the MS might there be pul- 
sating variables of t he Cephei or Sl owly Pulsating B 
stars types (see, e.g.. lDiago et al.ll2008l ). A sample of the 
brightest, intermediate, and the faintest light-curves is 
shown in Fig. [14J 

10. PREVIOUSLY KNOWN VARIABLES 

The variable stars in our ACS field in IC 1613 that 
were already known prior to this work are presented in 
Table [3 We give their index and period as measured 
from our data, as well as their identifications and peri- 
ods in previous catalogs; most of them were discovered 
during the wide-field survey by the OGLE collaboration. 
Because of the relatively shallow magnitude limits of the 
other surveys, all are luminous Cepheids with relatively 
large amplitude. 

Two of them, V1 24 and V172, were already included 
in Baade's catalog (Sandagel ll97lD as Vll and V35, re- 
spectively. These are the two brightest Cepheids of our 
sample, although the very long period of the former pre- 
vented us from measuring its intensity-averaged magni- 
tude and period. Another candidat e variab l e tha t was 
flagged as a possible Cepheid in iSandagel (1 1 9 7 If) and 
j udged to be an irregular variable in lCarlson fc Sandagd 
(1990), V41 is constant within ^0.02 and 0.05 mag in 
our F475W and F814W data and was not included in 
our catalog. 

Two other Cepheids appeared in iMantegazza et all 
(2001) with periods similar to the ones found here within 
0.01 day, while three LPV/Irr variables of their cata- 
log appear constant in our data (V0762C, V1598C, and 
V1667C). However, all have very long periods (>70 days) 



TABLE 7 

Cross-Identification with Previously Known Cepheids 
in IC 1613 



ID 


Period 


I T ) r\ f t 


POGLE 


TT~} j. i» ^ 

otner 


f other 


V008 


1.69 


2124 


1.697 


V1478C 


1.701 


V022 


1.48 


2197 


1.459 






V057 


1.69 


2818 


1.661 






V077 


1.35 


2909 


1.309 






V086 




2389 


2.029 






V093 


1.75 






V1465C 


1.741 


V106 


2.13 


2117 


2.131 






V124 




1987 


25.862 


Vll 


25.7719 


V137 


~3. 


2751 


1.376 






V147 


2.0 


2342 


1.972 






V154 


2.7 


2760 


2.712 






V170 


1.33 


2771 


1.329 






V172 


3.0 


2240 


3.074 


V35 


3.073417 



a Vll fc V35 from TSandagel {TWuh , V1478C & V1465C from 
Mantcgazza ct al. (2001) 



and/or low amplitude (Avj//i=0.20) and their variation 
could have easily been missed in our data due to the 
short observing run. 

The r emaining variables were previously observed by 
OGLE (jUdalski et al.1 120011 ). For most of them, the pe- 
riods determined from our data agree well with theirs. 
However, the 1.376 day period given by OGLE for V137 
cannot phase our data. Our sampling over 3 consecutive 
days clearly covers only about one cycle, so the period 
must be very close to P=3 days. 

By comparing the number of Cepheids observed by 
us and the OGLE team in the same field, we can ob- 
tain a rough estimate of the completeness of their sur- 
vey. We found 49 Cepheids over the area of the ACS 
camera, while their sensitivity limited them to 12 detec- 
tions. Assuming that the density of Cepheids is uniform 
in the central region of IC 1613, this means that about 
500 Cepheids should be present in the 14'xl4' field ob- 
served by the OGLE collaboration. Given the extent of 
the galaxy, this value is roughly in agreement with the 
total number estimated in § [4] 

11. DISTANCE 

In iPaper 1 we used the photometric and pulsational 
properties of the RR Lyrae stars to measure the distances 
to the dSphs Cetus and Tucana, and showed that the val- 
ues we obtained are in excellent agreement with the dis- 
tances obtained previously with independent methods. 
We therefore obtain here a RR Lyrae-based distance for 
IC 1613 in the same way. Additionally, given the signif- 
icant number of Cepheids we found in IC 1613 and the 
high quality of their light-curves, together with their firm 
classification as fundamental, first-, or second-overtone 
pulsators, it is possible to calculate accurate distances 
from their properties as well. 

In this section, we use several methods adopted in the 
literature to calculate the distance based on the pho- 
tometric and pulsational properties of the Cepheids and 
RR Lyrae stars. In all cases, the intensity-averaged mean 
magnitudes in the Johnson bands are used. 

11.1. Cepheids 

11.1.1. Period- Wesenheit relation for Fundamental and 
First-Overtone Cepheids 
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The main distance indicator for star forming galaxies 
in the nearby universe is the relation between the period 
and the luminosity of their Cepheids, even though the 
slope and zero-point of this rel ation might be s l ightly 
dependent on m e tallicity (e.g., IKennicutt et al.l 119981 : 
ISakai et al.l 120041: ISandage fc Tammannl I2008D . How- 
ever, fro m the analysis of nonlin ear convective pulsations 
models, iFiorentino et al.l (|2007l ) showed that the I-band 
period-luminosity is not expected to show metallicity de- 
pendence at metallicities lower than that of the LMC 
(Z=0.008), and that the effect of metal-content on the 
P-My relation is negligible below Z=0.004. Given the 
low metallicity of the young stars in IC 1613 (Z^0.003), 
we will therefore assume that it is safe to estimate its dis- 
tance by comparison with the properties of the Cepheids 
of the Magellanic Clouds. 

Here we chose to derive a new P-Wj (see § relation 
based on the Cepheids of both MCs instead of using the 
relations available in the literature for the following rea- 
sons: i) first, the domain of validity for the PL relations 
of the lit erature starts at about P=2.5 days (log P = 0.4; 
see e.g.. lUdalski et al.lll999cT: IFouaue et al.ll2007fl . while 
all our well-measured Cepheids have periods shorter than 
this value; and ii) the period distribution of the Cepheids 
in the LMC and SMC are very different from each other, 
with the peaks in the distribution at P^3.2 and 2.1 days 
for the fundamental and first-overtone Cepheids in the 
LMC, and at P~1.6 and 1.3 days, respectively, in the 
SMC. While the period distribution of the Cepheids in 
IC 1613 (with the peaks at P~1.6 and 1.1 days) are very 
similar to those in the SMC, the distances to other galax- 
ies are usually provided relative to a given LMC dis- 
tance. In addition, the PL relation of the SMC presents 
a rather large dispersion due to the inclination of the 
galaxy with respect to the line of sight. We thus com- 
bined the Cepheids of both Magellanic Clou ds assuming 
A(m-M) =0.51±0.03 (jUdalski et al.lll999d ) to improve 
the period coverage. This also has the advantage of can- 
celling any small difference in slope that might be present 
between the PL relations of the LMC and SMC. 

In addition, most of the Cepheids we observed in 
IC 1613 are located at the short-period end of the dis- 
tribution, so any non-linearity in the PL relation would 
lead to a different dist a nce h ad the Cepheids had longer 
periods. iNgeow et al.l f|2009T ) showed that the PL rela- 
tion in the Wesenheit magnitude Wj is linear over the 
period range covered by the LMC Cepheids, while it is 
not in the V&I bands. Finally, the use of Wi also limits 
the effect of interstellar reddening and thus reduces the 
scatter in the relation. 

We calculate the P-Wi relation from linear regression 
fits to the combined MC Ceph eids from OGLE-II ob- 
servations (jUdalski et al.lfl999bTldh over the whole range 
of periods after rejecting the outliers through sigma- 
clipping (2.5-ct, 5 iterations). Assuming an LMC dis tance 
of (m-M) =18.515±0.085 (jClementini et al.|[200l . we 
found: 

W/ F =-3.435(±0.007) log P F - 2.540(±0.006), and 
W/^ o =-3.544(±0.013) log P FO ~ 3.067(±0.007) 
for the fundamental and first-overtone Cepheids, respec- 
tively, with a standard deviation of 0.096 in both cases. 

The slope of our fundamental P—Wi relation is 
marginally different from the slope fou nd by other 
authors (e.g., W/ iF /logPF=3.313±0.008, INgeow et all 
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Fig. 15. — Distribution of RR Lyrae in the (V) -log P plane, 
where the predicted edges of the instability strip of lCaputo et alj 
(2000) have been overplotted (see text for details). Filled circles 
and triangles represent RRaft and RRc, respectively. 

12001 Wr F /logP P =3.320±0.011. IFouaue et~aTI l200l 
based on the OGLE Cepheids of the LMC only, 
and may be an indication that this relation deviates 
from linearity at periods shorter than about 0.3 days. 
On the other hand, it is intermediate between the 
slopes found for the LMC and t he values obtained 
for the Cepheids of the MW (e.g.. IFouaue et al.l [20071 
W FF /\ogP F =3.477±0.074). 

We applied the P—Wj relations derived above to 
the fundamental and first-overtone Cepheids of IC 1613 
shown in Fig. |H] (excluding the two Cepheids located at 
the midpoint between the dash-dotted and dashed lines, 
V008 and V093), and found (m-M) =24.50±0.11 and 
24.47±0.12, respectively. The good agreement between 
these two values indicates that first-overtone Cepheids 
can be safely used as distance indicators in the case where 
the pulsating mode of the Cepheids can be unambigu- 
ously identified. 

We also calculated the distance to IC 1613 using the 
PL relation f or fun damental-mode Cepheids given by 
INgeow et al.1 (|2009D to check if the small differences 
in slope and zero-point described above would have 
a significant effect, but found a very similar value 
((m— M)o=24.55±0.12) within the uncertainties. 

11.1.2. Second-Overtone Cepheids and the 
Period-Luminosity- Color relation 

From t heir theoretical wo rk on second-overtone (SO) 
Cepheids, IBono et ail ((2001) noted that these variables 
follow a period-luminosity-color relation of the form: 



-3.961 - 3.905/O5P + 3.250(y - I). 



Given the low amplitude of the luminosity variations and 
the very small temperature width of the SO IS, the stan- 
dard deviation of this relation is only 0.004. It is thus 
possible to use it to determine distances with good ac- 
curacy. From the SO Cepheids observed in the SMC by 
lUdalski et all ()1999al ) and assuming E(B— V)=0.054, the 
authors found a distance of (m— M) ,smc=19.11±0.08 to 
the SMC through this relation. 

Applying this relation to the two SO Cepheids de- 
scribed in § 16.21 we find a distance modulus of 
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TABLE 8 

Summary of distance moduli to IC 1613. 
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F PL 
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F PL 
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F PL 
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F PL 
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Wj 
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F PL 
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F PL 
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F PL 
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18.50 


.1 
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F PL 
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- 


24.43±0.08 
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F PL 
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— 


1 
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1 
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0.025 
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F PL 
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F PL 
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F PL 
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0.09 
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F PL 
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24.50±0.12 


0.024 


18.54 
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F PL 
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24.43±0.05 


0.025 




I 


TRGB 
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24.29±0.07 


0.08 


18.50 


3.6/im 


Cep: F PL 


14 


24.28±0.07 


0.08 


18.50 


4.5/im 


Cep: F PL 
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24.0 



24.5 



Cep: F PL (B) 


Sandage 1971 


Cep: F PL (B) 


de Vaucouleurs 1978 


Cep: F PL (H) 


McAlary+1984 


Cep: F PL (BVRIH) 


Freedman 1 988 


Cep: F PL (W,) 


Freedman 1 988 


RRL: M v -[Fe/H] 


Saha+1992 


Cep: F PL (V) 


Macri+2001 


Cep: F PL (1) 


Macri+2001 


Cep: F PL (J) 


Macri+2001 


Cep: F PL (H) 


Macri+2001 


TRGB 


Dolphin+2001 


Red Clump 


Dolphin+2001 


RRL: M u -[Fe/H] 


Dolphin+2001 


Cep: F PL (V) 


Udalski+2001 


Cep: F PL (1) 


Udalski+2001 


Cep: F PL (W,) 


Udalski+2001 


TRGB 


Tikhonov+2002 


Cep: F PL (J) 


Pietrzynski+2006 


Cep: F PL (K) 


Pietrzynski+2006 


Cep: F PL (BVRI) 


Antonello+2006 


TRGB 


Rizzi+2007 


Cep: F PL (3.6 (im) 


Freedman+2009 


Cep: F PL (4.5 Jim) 


Freedman+2009 


Cep: F PL (W,) 


this work 


Cep: FO PL (W,) 


this work 


Cep: SO PLC 


this work 


RRL: M v -[Fe/H] 


this work 


RRL: PLM 

, 1 , , , , 


this work 


25.0 25.5 


26.0 


(m-M) 





Re fere nces. — (1) This w ork; (2) ISaada.geH1971l : ( 3) Ide Vaucouleur j 
11973 (4) IMcAlarv et al.l 119841: (5) IFreedmarill98a t '6l ISaha et all 1 199? 
(7) Maori et all 120011: fsTlDorphin et al.l I200ll: (9) lUdalski et all booil : 
(10) Tikhonov & Galazutdinova 2002; (11) Pictrzyriski ct al. 2006; (12) 
lAntonello et al.1120061 ; (13) IRizzi et al.N2007l ; (14) IFreedman et al.H200a 

(m-M) =24.54±0.11. Note, however, that the distance 
to the SMC used in the previous section to obtain the 
period- Wj relationship was 19.03, assuming an LMC dis- 
tance of 18.515±0.085 and a distance moduli difference of 
0.51±0.03. Taking into account this offset in zero-point, 
we find that the true distance modulus to IC 1613 from 
SO Cepheids is 24.46±0.11, in better agreement with the 
mean value determined below. 

11.2. RR Lyrae stars 

Similar to what we did in iPaper 1 we also calculate 
the distance to IC 1613 using two methods based on the 
properties of the RR Lyrae stars, namely the luminosity- 
metallicity relation, which arises from the knowledge that 
the intrinsic luminosity of HB stars mainly depends on 
their metallicity, and the period-luminosity-metallicity 
(PLM) relation, based on the theoretical location of the 
instability strip in the period-luminosity plane. 

The luminosity- metallicity relation we used in IPaper II 
has the form: 

M v = 0.866(±0.085) + 0.214(±0.047)[Fe/ii]. 

To calculate the mean magnitude of the RR Lyrae stars, 
we used only the stars for which we could determine ac- 
curate intensity-averaged mean magnitudes — RRa6 and 
RRc — and obtained (V r )=24.99±0.01. Assuming a mean 
metallicity of Z=0.0005 ±0.0002 (i.e., [Fe/H]=-1 6±0.2 
assuming Z =O.O198 ((G revesse fe Noels! 11991 and 
fFe/H l=log Z+1.70 - log(0.638 f + 0.362) (|Salaris et al.1 
1993) with the a— enhancement factor f set to zero) 
for the old population from the results of the SFH dc- 



FlG. 16. — Summary of the distance moduli from the present 
work and the literature. The reddened moduli as given in the ref- 
erenced articles are shown as pluses, while the moduli corrected for 
a common reddening of E(B— V)=0.025 and LMC distance mod- 
ulus (m— M)o=18.515±0.085 are shown as filled circles. The ver- 
tical lines indicate the mean and st andard deviation ( ^=0.071) of 
these measurements (excluding the ISaha e t al. (1992) outlier) at 
(m-M) =24.400±0.014 (statistical). 

rived in Skillman et al. (2010, in preparation), we find 
a luminosity for the HB of M y =0.52±0.12. The un- 
certainty was quantified through Monte-Carlo simula- 
tions. This gives a distance modulus of 24.47±0.12, or 
(m— M)o=24.39±0.12 after correcting for reddening. 

The second method we used to determine the dis- 
tance modulus consists of matching the PLM relation 
at the first-overtone blue edge (FOBE) of the IS from 
iCaputo et al.l (|2000l their eq. 3) to the observed RRc. 
Basically, the theoretical limits of the IS are shifted in 
magnitude until the FOBE coincides with the observed 
distribution of first-overtone RR Lyrae stars. Figure [T5l 
shows the position of the IS {dashed lines) overplotted 
on the distribution of observed RR Lyrae stars in the 
{V)-\og P plane. 

Adopting the mean metallicity given above 
(Z=0.0005), we derive a deredd ened distance mod- 
ulus of (m-M) =24.36, for which ICaputo et aLI (pOOOf ) 
give an intrinsic dispersion of ov=0.07 mag due to 
uncertainties associated with the various ingredients of 
the model. Combined with the errors in metallicity, 
mean magnitude, and period, we estimate a total 
uncertainty in this measure of the distance of ~0.1. 

11.3. Results 

Since the distances derived in § 111.11 and 111.21 are con- 
sistent with each other and have similar uncertainties, 
we simply averaged them to obtain our best distance es- 
timate and obtained a dereddened distance modulus of 
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24.44 (er=0.054), corresponding to a distance of 770 kpc. 

A summary of the previous distance determination 
from the literature is given in Tabled The columns give 
the reddened distance modulus, the assumed or calcu- 
lated reddening and LMC distance modulus, the observa- 
tion band, and the method used to estimate the distance. 
We compare these values in Fig. [THl where the reddened 
moduli as given in the referenced articles are shown 
as pluses. To make this comparison more meaningful, 
we also show the distance moduli corrected for a com- 
mon reddening of E(B-V)=0.025 (jSchlegel et all 11998ft 
and LMC distance m odulus (m-M) =18.515±0.085 
(jClementini et all [2003V The red dening correction was 
done following the extinction law of lCardelli et al.l (|1989l ) 
with Ry=3.1. 

The figure shows that all the measurements, indepen- 
dent of the method, are in excellent agreement with each 
other once s et on the same sca le. From these values (ex- 
cluding the iSaha et""aH (|1992h outlier), we find a mean 
dereddened distance of (m-M) =24.400±0.014, where 
the uncertainty only includes the standard deviation of 
the different measurements. The vertical lines in Fig. 1161 
indicate the mean and standard deviation (<t=0.071) 
of these values. Therefore, we believe that the largest 
source of uncertainty on the distance to IC 1613 is sys- 
tematic rather than statistical, and mostly lies in the 
calibration of the distance to the LMC. 

12. DISCUSSION AND CONCLUSIONS 

We have presented the results of a new search for vari- 
able stars in IC 1613 based on high quality HST images. 
In the ACS field, we found 259 candidate variables, in- 
cluding 90 RR Lyrae stars, 49 Cepheids, and 38 eclipsing 
binaries, as well as nine RR Lyrae stars and two addi- 
tional Cepheids in the parallel WFPC2 field. Only thir- 
teen of these variables were known prior to this study, 
all of them Cepheids. 

We find that the mean periods of the RRafe and RRc 
stars, as well as the fraction of overtone pulsators, place 
this galaxy in the intermediate regime between the Oost- 
erhoff types, as was already observed for the vast major- 
ity of Local Group dwarf galaxies. 

From the comparison with Magellanic Clouds Cepheids 
through the period-luminosity diagram, we find that all 
our Cepheids are bona fide classical Cepheids, while the 
Cepheids in the dSph Cetus and Tucana discussed in Pa- 
per I were classified as anomalous Cepheids. The lack of 
classical Cepheids in the latter dwarfs is not surprising 
given the absence of signifi cant star formation more re- 
cent than about 8-10 Gyr (|Monelli et al.ll2010D . On the 
other hand, the reason for the apparent lack of anomalous 
Cepheids in IC 1613 is not clear. T he star formation his- 
tories calculated from deep CMDs (jSkillman et al.ll2003L 
2010, in preparation) indicate that it formed stars fairly 
constantly over the past 13 Gyr, so unless the formation 
of anomalous Cepheids is extremely sensitive to the metal 
content, anomalous Cepheids should also be present in 



I C1613. 

IGallart et"aTl (|2004h showed that both types of 
Cepheids were present in comparable numbers in the 
transition type galaxy Phoenix, and possibly also in Leo 
A and Sextans A, all of them having very low metal- 
licity (Z<0.0008). These authors suggest that the only 
requirements for the presence of both types of Cepheids 
in a galaxy are low enough metallicity (Z~0.0004) and 
star formation activity at all ages. At the other end 
of the dwarf galaxy metallicity spectru m, 83 anomalous 
Cephe ids were discovered in the LMC (fSoszvhsk i et al.l 
2008b), compared to the 3361 classical Cepheids tha t 
are present in the same catalog f Soszvhs ki et al.ll2008al) . 
This implies about 2.4 anomalous Cepheid per 100 classi- 
cal Cepheid at the metallicity of the LMC (Z~0.008). If 
this fraction is representative for low metallicity galaxies, 
about one anomalous Cepheid could be e xpected in our 
ACS field-of-view of IC 1613. However, Dolphin et al. 
( 2002) showed how low metallicity galaxies contain many 
more short-period classical Cepheids than higher metal- 
licity galaxies at a given star formation rate due to mor- 
phology of the blue loops: at low metallicity, the blue 
loops extend further to the blue and thus cross the in- 
stability strip at fainter magnitudes. As a result of the 
stellar initial mass function, the fainter Cepheids are also 
more numerous. Therefore, the apparent lack of anoma- 
lous Cepheids in our field might simply be due to small 
number statistics. Given the low density of anomalous 
Cepheids expected in IC 1613, more data covering a large 
field-of-view are needed to check if anomalous Cepheids 
are actually present in this galaxy. 

Finally, we used the properties of the RR Lyrae stars 
and Cepheids to estimate the distance to IC 1613, and 
found excellent agreement with the values previously de- 
termined. Combining all the measurements after correc- 
tion for a common reddening and reference LMC dis- 
tance, we find a true distance modulus to this galaxy 
of (m-M) =24.400±0.014 (statistical), corresponding to 
760 kpc. 
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APPENDIX 

A. COMMENTS ON INDIVIDUAL VARIABLES 

The following comments arc based on the careful inspection of the stacked images, of the light curves, and/or 
peculiar properties exhibited on one or more of the plots presented in this work. 
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V001 
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Blend. 


V004 
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Blend. 


V007 
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Close to a bright star. 
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Blend. 
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Possible blend. 
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Close to the wing of a very bright star. 
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Blend. 


V036 


— 


Blend and close to a bright star. 
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Blend. Very low amplitude RRafc. 
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Blend. 
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Blend. 
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Close to a very bright star. 
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Blend. 
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Blend. 
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Blend with a background galaxy and bright stars. 
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Blend. 


V077 
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Blend. 
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In bad column. 


V085 
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Blend. 
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Blend. 
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Blend. 
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Possible blend. 
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Blend. 
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Blend. 
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Located close to the edge of chip 1. Some points bad or missing 
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Blend. 
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Blend. 
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Blend. 
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Possible blend. 


V121 




Blend and located close to the edge of chip 1. Some points bad c 
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Blend. 


V125 
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Blend. 


V133 




Blend. 
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Blend. 


V144 
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Blend. 
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Possible blend. 


V148 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 


V171 




Blend. 


V175 




Possible blend. 


VI 76 




Blend. 


V177 




Blend. 


VI 78 




Blend. 


V180 




Blend. 



Faint Variables in IC 1613 



19 



VC202 - 


Blend. 


VC206 - 


Blend. 


VC208 - 


Saturated in F814W. 


VC213 - 


Possible blend. 


VC217 - 


Blend. 


VC223 - 


Blend. 


VC225 - 


Possible blend. 
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Blend. 


VC229 - 


Blend. 


VC231 - 


Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Blend. 
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Possible blend. 
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Possible blend. 
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Possible blend. 


VC252 - 


Blend. 
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Blend. 


VC257 - 


Possible blend. 


VC259 - 


Blend. 


VC265 - 


Blend. 


VC267 - 


Blend. 


VC268 - 


Blend. 


VC270 - 


Blend. 



B. FINDING CHARTS 

The finding charts for the whole sample of variable stars are presented in the electronic version of The Astrophysical 
Journal (Fig. IT?] . 
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